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ABSTRACT 

Planets acquire atmospheres from their parent circumstellar disks. We derive a general analytic 
expression for how the atmospheric mass grows with time t, as a function of the underlying core mass 
Mcore and nebular conditions, including the gas metallicity Z. Planets accrete as much gas as can 
cool: an atmosphere’s doubling time is given by its Kelvin-Helmholtz time. Dusty atmospheres behave 
differently from atmospheres made dust-free by grain growth and sedimentation. The gas-to-core mass 
ratio (GCR) of a dusty atmosphere scales as GCR ex where ft^ch oc 1/(1 — Z) (for 

Z not too close to 1) is the mean molecular weight at the innermost radiative-convective boundary. 
This scaling applies across all orbital distances and nebular conditions for dusty atmospheres; their 
radiative-convective boundaries, which regulate cooling, are not set by the external environment, but 
rather by the internal microphysics of dust sublimation, H2 dissociation, and the formation of H“. By 
contrast, dust-free atmospheres have their radiative boundaries at temperatures Trcb close to nebular 
temperatures Touti and grow faster at larger orbital distances where cooler temperatures, and by 
extension lower opacities, prevail. At 0.1 AU in a gas-poor nebula, GCR oc 

while beyond 1 AU in a gas-rich nebula, GCR oc We confirm our analytic 

scalings against detailed numerical models for objects ranging in mass from Mars (O.IM0) to the 
most extreme super-Earths (IO-2OM0), and explain why heating from planetesimal accretion cannot 
prevent the latter from undergoing runaway gas accretion. 


1. INTRODUCTION 

The Kepler mission has discovered that at least ^50% 
of Sun-like stars harbor “super-Ea rths”—here defined as 
planets having radii l-4i?0 (e.g., Fressin et al. 2013) ^ 
Their masses, measured by tr ansit timing variations 
(e.g., Hadde n fc Lithwick 20 14|) a nd Doppler radial ve- 
locities (e.g., Weiss &: Marcy||20l4 1, imply bulk densities 
that are typically < 3 g cm These densities are too 
low to be compatible wi th a pure rock composition. The 
consensus view (see also Rogers|[2015 ) is that the masses 
of super-Earths are doininated by their solid cores — 
of mass Mcore — 2-2OM0 and radius i?core — l-2i?0 
— while their total radii can be more than doubled by 
voluminous, hydrogen-rich atmospheres. Interior mod¬ 
els su ggest gas-to-core mass r atios (GCRs) up to ^10% 
(e.g., Lopez fc Fortney 2014) and more typically ^1% 

( Wolfgang &: Lopez||2U15 ). 

Unlike the Earth’s atmosphere (GCR ^ 10“®), the at¬ 
mospheres of extrasolar super-Earths are likel y too mas¬ 


sive to have been outgassed from rock (e.g., Rogers & 
Seager 2010). More plausibly, super-Earth atmospheres 


originated as the envelopes of gas giants like Jupiter 


Ikoma & Hori 

2012 Boden- 

Inamdar & Sch 

ichting 
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hnd that super-Earth cores can accrete atmospheres hav- 
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ing GCRs of ~1-10% before the gas disk dissipates on 
Myr timescales. Even higher GCRs can be obtained 
under a variety of conditions (Lee et al. 2014). These 
higher values may be required because once the parent 
nebula disperses and planets are laid bare, atmospheric 


loss driven by stellar irr adiation (e.g., Lopez & Fortney 
2013 Owen fc Wu|2013) and by the young planet's heat 


of formation (|Owen & Wu 2015) can pare GCRs down 
b y factors of severa l or more. 


Lee et al. (2014), hereafter paper I, computed nebu¬ 
lar accretion histories for a range of core masses, disk 
temperatures and densities, and atmospheric metallici- 
ties and dust contents. Our aim here, in paper H, is 
to provide an analytic understanding of their numeri¬ 
cal results. Benefiting from hindsight, we will reduce 
our model to a few essential elements and obtain sim¬ 
ple power-law scalings between GCR, time t, core mass 
Mcore, and metallicity Z. These scalings will be derived 
against a variety of backdrops: gas-rich vs. gas-poor neb¬ 
ulae; dusty vs. dust-free atmospheres; close-in vs. far-out 
orbital distances. 

Before we present these scaling relations (Section]^, 
we revisit the fundamental assumption underlying them: 
that the nascent atmospheres have no power source but 
passively cool by radiating into their nebular wombs. 
One source of power that we ignore — but that is com¬ 
monly invoked in the literature — is the accretion of 
planetesimals. Paper I provided reasons why planetesi¬ 
mal accretion could plausibly be dropped when consid¬ 
ering the origin of super-Earth atmospheres; in Section 
o below, we flesh these arguments out more fully and 
quantitatively. Readers interested in the main results 
of this paper can skip ahead to Section which derives 
the GCR (t, Mcore, Z) scalings for passively cooling atmo- 























































2 


spheres, and to Section|^which contains a recapitulation 
with commentary. 


1.1. Planetesimal Accretion 

Paper I articulated one of the main puzzles posed by 
super-Earths: how, despite their large core masses, they 
avoided being transformed into Jupiter-mass giants by 
runaway gas accretion, and instead had their GCRs sta¬ 
bilized at values of several percent. 

One way to stop an atmosphere from growing is to sup¬ 
ply it with sufficient heat — enough to bafance cooling 
and arrest secular contraction. The energy released by 
the accretion of planetesimals is a candidate heat source. 
We can estimate the required rates of mass delivery M 
by equating the accretion luminosity 


RrorR 


( 1 ) 


to the cooling luminosity Tcooi, where the latter is com¬ 
puted from the models presented iripaper I. Here G is the 
gravitational constant. Equation (IH assumes that plan¬ 
etesimals are large enough to penetrate the atmosphere 
and release their kinetic energy at the core surface. 

Figure llj displays the planetesimal accretion rates M so 
estimated^ using the dusty atmosphere models from pa¬ 
per I (dust-free atmospheres will be considered shortly). 
Every curve exhibits a minimum in M with GCR; this 
minimum corresponds to the minimum in Lcooi with time 
that appears in all passively cooling a nd growing atmo¬ 
spheres (see also Piso & Youdin 2014) and that we used 
in paper I to mark the onset oi runaway gas accretion. 
To the left of the minima, the increasing atmospheric 
density with increasing GCR renders the envelope more 
opaque and causes Tcooi to drop (see Section 3.1 of pa¬ 
per I and Section!^ of this paper). To the right of the 
minima, at GCR ^0.5, the self-gravity of the gas enve¬ 
lope becomes significant, and larger Lcooi is required to 
balance stronger gravity. 

The curves in Figure represent the loci of thermal 
equilibrium; atmospheres are stabilized at a given GCR 
when solids rain down at the corresponding M. We now 
consider the stability of these equilibria. To the left of the 
minima, at lower GCR, equilibria are stable; for example, 
a perturbation to higher GCR decreases Tcooi below Lace 
(which is presumed fixed), resulting in a net heating that 
expands the atmosphere and lowers the GCR back down 
to its equilibrium value. To the right of the minima, 
at higher GCR, there are no stable equilibria because 
Lcooi increases with both GCR and Mcore- As Mcore 
increases from planetesimal accretion, Lcooi rises rapidly 
(see Section 3.2.3 of paper I and Sectionj^of this paper), 
outpacing Lace and triggering runaway. 

The values of M to the left of the minima in Figure]^ 
represent possible solutions to the puzzle of how super- 
Earths avoided runaway gas accretion. But we find these 
solutions unsatisfactory for a couple of reasons. The first 
is that planetesimaf accretion rates must be fine-tuned to 
the values plotted. Why, for example, should 5Mf^ cores 
accrete O.1M0/Myr in solids to have their GCRs stabi¬ 
lized at ~0.1? Planetesimal accretion rates are influenced 


by a host of factors (e.g., Goldreich et al. 2004), and most 
estimates lead to rates orders of magnitude higher than 



Figure 1. Planetesimal accretion rates required to keep dusty 
envelopes in thermal equilibrium at a given gas-to-core mass ratio 
GCR. Mass infall rates M are calculated by equating the accretion 
luminosity Lace to the cooling luminosity Lcooi 5 where the latter 
is computed from our numerical models of dusty, passively cooling 
atmospheres in a minimum-mass extrasolar nebula (MMEN) at 0.1 
AU (paper I; for reference, the gas surface density S of the MMEN 
is several times larger than that of the traditional minimum-mass 
solar nebula; the difference is immaterial for all the results pre¬ 
sented in this paper). Circles mark minima in M which correspond 
to minima in Lcooi I these minima, GCR ~ 0.5 and envelopes are 
on the brink of runaway gas accretion. To the left of these minima, 
at smaller GCRs, are stable equilibria (thick lines) and to the right 
are unstable equilibria (thin lines). To stabilize GCRs at values 
< 0.5 requires that planetesimal accretion rates be fine-tuned to 
the values plotted. In the case of cores of mass IOMq, this fine- 
tuning would still be unable to prevent runaway: the required M’s, 
of order Myr“^, would double the core mass within a disk 

lifetime of ~10 Myr and push the cores over to runaway (see also 
Figure §. 


the ones ploH ed in Figure IT] (see, e.g.. Appendix A of 
Rafikov||2006 ). The much hi^i er accretion rates are nat- 


ural consequences of the high disk surface densities and 
short orbital times characterizing the small orbital dis¬ 
tances where Kepler super-Earths are found, and lead us 
to surmise that planetesimals are fully incorporated into 
planets on timescales much shorter than the 1-10 Myrs 
required for atmospheres to cool and growj^ 

The second reason we do not find planetesimal accre¬ 
tion attractive is that it bites off more than it can chew: 
Figure n implies that a IOMq core has to accrete plan- 
etesim^ at such a high rate that it doubles in mass 
within the disk lifetime of tdisk ~ 10 Myr. The 2OM0 
core that results is certain to undergo runaway gas ac¬ 
cretion in a disk that is gas-rich (i.e., with a gas surface 
density E comparable to that of the minimum-mass solar 
nebula). 

We can illustrate with another calculation how plan¬ 
etesimal accretion leads to problems of fine-tuning and 
runaway gas accretion. An atmosphere does not grow 
if Aacc > Lcooi- In particular, runaway is avoided if 


^ Our discussion here pertains to “planetesimals”, not to the 
larger protocores which are thought to assemble into super-Earth 
cores by “major mergers”. Such mergers take place over a wide 
range of timescales that can overlap or even exceed gas disk life- 
times. depending se nsitively on the disk’s solid surface density 
||Dawson et al.|2015||. Energy release from the last major merger 
could, m principle, provide a significan t source of heat to stop gas 
accretion pnamdar &: Schlichting|2015|l, but only if the core trans¬ 
ports its internal energy outward on a timescale comparable to the 
atmospheric cooling time of 1-10 Myr. The actual energy trans¬ 
port timescale of the core is highly uncertain, depending on the 
unknown viscosity (see paper I, section 3.1.2). 
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the scaling relation between Lcooi and Mcoie as derived 
in Section below.) Figure plots Facc/icooi|GCR= 0.5 
against M. Cores of initial mass M^ore = can avoid 
runaway — but only if their atmospheres are dusty, and 
if M is tuned to a narrow range. This window closes 
completely for cores of initial mass > 7M^: such cores 
will inevitably undergo runaway within 10 Myr in gas- 
rich disks, regardless of the magnitude of planetesimal 
accretion or their proximity to their central stars. These 
conclusions are only amplihed for cores with dust-free 
atmospheres. We will elaborate in Section on these 
trends with nebular environment and atmosjmeric com¬ 
position. The point here is that planetesimal accretion 
does not generically prevent runaway. Avoiding runaway 
for even low-mass cores with Mf-ore < SMq requires del¬ 
icate adjustment of M which seems difficult to achieve 
under general circumstances. We are motivated there¬ 
fore to study gas accretion histories that omit planetes¬ 
imal accretion — this is the main subject of this paper, 
to which we now turn. 

2. SCALING RELATIONS FOR OCR (Mcore, b Z) 

We derive general scaling relations for how the gas-to- 
core mass ratio GCR varies with core mass Mcore, time 
t, and metallicity Z, in the absence of external heating. 
The derivation is semi-analytic in that a few parameters 
will be calibrated using our numerical models (see paper 
I for details on how we build our numerical models). 

Accretion is mediated by cooling: upon radiating away 
its energy, a planet’s atmosphere contracts, allowing 
nebular gas to refill the Hill sphere. The system self- 
regulates so that whatever atmosphere of mass Mgas has 
been accreted has a cooling time equal to the time that 
has elapsed: 


Mgas _ , _ , _ \77'\ 

71 V i ^COOl j 

Mgas ^cool 


( 3 ) 


Figure 2. Preventing runaway with heating from infalling plan- 
etesimals is possible only for specific ranges of mass infall rates M 
for sufficiently low-mass cores. Accounting for the growth of the 
core mass from planetesimal accretion (curves are labeled by initial 
core masses), we plot the ratio of Lace to Lcooi (GCR = 0.5) after 
a gas lifetime of tdisk = 10 Myr in a gas-rich (MMEN) disk, con¬ 
sidering both dusty and dust-free atmospheres from 0.1 to 5 AU 
(see annotation in each panel). If Lacc/Lcooi ^ 1 (unshaded), then 
runaway is successfully avoided (see equation and surrounding 
discussion). All four panels show that the ranges of M and initial 
core mass required to avoid gas giant formation are extremely lim¬ 
ited — if they exist at all. If M is too low. Lace can never balance 
Lcooh if M is too high, cores rapidly grow in mass and become 
giants within the disk lifetim e. The lowest M’s estimated from first 
principles by |Rafikov| | |200^ see his equation Al) are indicated by 
red lines; these M’s readily push cores over to runaway. 

Lacc > Lcooi when GCR ~ 0.5, the critical GCR above 
which envelope self-gravity starts becoming signiheant. 
We therefore evaluate the ratio L^cc/Lcooi when GCR = 
0.5, presuming the core has accreted planetesimals at a 
constant rate for the preceding time tdisk^ 


.bcool 


GCR=0.5 


G{McoTe + Mtdlslc)M 
L^coceLcool 


( 2 ) 


where Tcooi is evaluated using our atmospheric models for 
the final core mass Mcoie + Mtdisk- (This evaluation uses 


where E is the atmosphere’s total energy and Tcooi is its 
luminosity. Statement (|^ is perhaps more easily under¬ 
stood by considering the inverse cases t <C tcooi (GCR 
is overestimated because not enough time has elapsed to 
accrete such a thick atmosphere) and t tcooi (GCR is 
underestimated because there is plenty of time for the 
atmosphere to continue cooling and growing). 

The relevant cooling time is that of the innermost con¬ 
vective zone which contains most of the atmosphere’s 
mass and energy. To estimate E, we use the fact that 
in hydrostatic equilibrium, an atmospheric mass A/gas = 
GCR X Mcoie has a total energy of order its gravitational 
potential energy: 


\E\ 


GMcoreAfgas 

R 


( 4 ) 


where G is the gravitational constant. What R should we 
choose: the core radius i?core or the radiative-convective 
boundary i?rcb"13 The answer depends on how steep the 
density profile is. For fixed adiabatic index 7 , the density 


® The outer radius of our numerical models is either the Bondi 
or Hill radius, whichever is smaller. Neither of these radii enters 
into our analytic theory, since not much mass is situated near the 
outer boundary. 
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profile in the isentropic convective zone follows 


P — Prcb 


1 + Vad ■ 


GM„ 


"-rcb 


1 


1 


-^rcb 




( 5 ) 


where prcb is the density at the radiative-convective 
boundary (rcb), Vad = (7 — l)/7 is the adiabatic gra¬ 
dient, = feTrcb/z^rcbTOH, T^ch and /Ltrcb are the tem¬ 
perature and mean molecular weight evaluated at the 
rcb, k is Boltzmann’s constant, and toh is the atomic 
mass of hydrogen. Now because 1/r > l/i?rcb and 
GMcore/c^cb^ > GMcore/c^^b^rcb ~ 1 (the last equal¬ 
ity follows from hydrostatic equilibrium), equation ([^ 
can be approximated as 


P ~ Prcb V 


Rh. 


rcb 


ad 


l/(7-l) 


( 6 ) 


where i?b,rcb = GMcore/c^cb- Equation (|^ implies that 
if 7 < 4/3, then the atmosphere’s mass is concentrated 
near i?core rather than near i?rcb- We hnd that the con¬ 
vective zones of all our numerical models are indeed char¬ 
acterized by 7 < 4/3: the adiabatic gradient drops at 
temperatures exceeding 2500 K as energy is spent dis¬ 
sociating H 2 rather than heatiM the gas. Therefore we 
choose R = i?core in equation Hn 

1^1 ^ GM/„„ X GCR 


Rr 




1/3 


mUI X GCR 
JeMUI X GCR 


( 7 ) 


where Je = G(47rpb/3)^/^ and pb is the bulk density of 
the core (assumed constant for this paper; we neglect the 
small variation of co re d ensity with core ma ss; see, e.g., 
Valencia et al.||2006 and Fortney et al.||2007 1. 

We now examine Tcoob 'i'he rcb controls the rate at 
which the innermost convective zone cools. Very little 
luminosity is generated above the rcb (as verified in paper 
I, section 3.1.2), so we evaluate Tcooi at the rcb: 

_ 647rG(l -I- GCR)McoreO’T)),bMrcbTOHVad 

-^COOl — 0 7 V^/ 

Oftjprcb^rcb 

where a is the Stefan-Boltzmann constant and Urch is the 
opacity at the rcb. We parameterize the latter as 

A«rcb = ^^oiPrcb/po^iTrcb/Tof{Z/Zo)^ (9) 

where the various constants depend on microphysics 
which vary from case to case (details to be given in the 
subsections below). 

We relate prcb to GCR as follows. The total atmo¬ 
spheric mass in the inner convective zone is 


Mgas = 47r 


^p{r)dr 


- 47rprcb(Vadi?b.rcb)^/^'^"'^i?L”ry^^”'^ (10) 

^ This is contrary to equation (32) in paper I which mistakenly 
assumes the atmosphere’s mass is concentrated near -Rrcb instead of 
near Rcore- The correction lengthens the runaway time estimated 
in that equation by a factor of 10, bringing it into closer agreement 
with the numerical result cited there. 


where we substituted Then 


GCR X Mr, 


Prch 


( 11 ) 


47r(VadRb,rcb)l/('>'-l)-RL”ry^"'”’^' 
Substituting (|^ and (11) into ([^: 

647rGi7mH f ZoY YrlYrcb^^d^ + GCR)Me, 


-^cool 


1 +a 

Prch 


3fcKoPo “T’o”'" V ^ . 

44+“7r2+“GamH ' r;eb^A‘rcbVad(l + GCR) 

Z 


SkKoPo “TJ 


— arji—p 


GCR4+" 


X M-“ (Vadi?b.rcb)^i?Sre"-^^^'^“^ 


( 12 ) 


We re-write (12) in terms of Mcore, Z,p.rcb, and GCR: 


where 


/l = 


.t^cool ^ I ^ I -^rcb GCR4^“ 

X (prcbVad)^+^Mcotl^^^ (13) 


44+“7r2+“a 

3«opF^ 


( 3 y4+a)[i-3^ 

dTTpb/ 


( 14 ) 


When all the hydrogen is molecular, the mean molec¬ 
ular weight pL depends on Z as: 


1 


p - 

V = 

Y = 


0.5V-h 0.25V-k 0.06Z 
1 - Z 


(15) 


1.4 

0.4(1 - Z) 
L4 


where X and V are the hydrogen and helium mass frac¬ 
tions, respectively. The prefactor of 0.06 for Z corre¬ 
sponds to the co ntribution from atomic metals using 
the abundances of|Grevesse &: N oels] (|1993D; thes e abun - 
dances are the ones adopted by ll'erguson et ai.| (|2005|), 
whose opacities we use. 

Collecting ([^ and (13) into ^ yields 


/zy GCR^+- j(i-/^) 
h \Zo) 1 + GCR 


X T 


-3-b/3-bi^ 


rcb 


(Mrcb V ad) 


_1_J 


(16) 


which we invert to arrive at our desired relation for GCR 
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as a function of f, Mcore, Z, and /ircb, valid for GCR < ij^ 
GCR = / 


(§ ] 


X (MrcbVad)^+^T"cb^ 


(17) 


We have introduced a dimensionless fudge factor / which 
we will normalize against our numerical models. The 
parameters that vary most from one scenario to another 
are the opacity constants in equation (|^, and the rcb 
variables T^cb and Vad- All these input constants will be 
drawn from our numerical solutions. 

The following subsections examine how equation 0 
plays out in various formation environments. We con¬ 
sider dusty vs. dust-free atmospheres in gas-poor vs. gas- 
rich nebulae at small orbital distances and large. The 
plausibility of these scenarios is not assessed; that exer¬ 
cise is deferred to a later study (paper III). In “dusty” 
m odels, we assume the I SM-like grain size distribution 
of Ferguson et al. (20051, and in dust-free models, we 
assume all metals to be in the gas phase. By “gas-rich” 
we mean a nebula whose gas surface density S equals 
that of the minimum-mass extrasolar nebula (MMEN; 
see eq uation 12 of paper I; for comparison, the Hayashi 
1981 nebula is ~7x less dense), and by “gas-poor” we 
mean a nebula whose gas content is 200 x smaller (one 
whose gas mass equals its solid mass). Nebular tempera¬ 
tures are taken from equation (13) of paper I (Tout = 
{1000,400,200} K at orbital distances a = {0.1,1,5} 
AU). 


2 . 1 . Dusty Atmospheres 

Dusty atmospheres are high opacity atmospheres and 
tend to be convective in their upper layers. They cease 
being dusty at depths below which temperatures are 
high enough for dust sublimation. The disappearance 
of grains causes the opacity k to drop by two orders of 
magnitude; the sudden transparency opens a radiative 
window at depth. This radiative zone appears univer¬ 
sally over all core masses and orbital distances as long 
as the upper layers are dusty. The base of the radiative 
zone — i.e., the innermost radiative-convective bound¬ 
ary — is located where H 2 dissociates and H“ appears 
with its strongly temperature-sensitive opacity: 

k{H-) ~ 3 X 10"^ cm^ g-i (" ■, ■_ ^ — 

\ 10 ^ g cm 

T y Y z 

2500 k; 





2.1.1. Dusty and Gas-Poor from 0.1-1 AU 
Substitu ting T^ch = 2500 K and the parameters from 


(18) into (17), and further restricting our attention to 
gas-poor nebulae for which the gas surface density E is 
1/200 that of the MMEN, we find 


GGR ~ 0.06 


1.2 

3.4 


t 


/2500K^ 


IMyry V '^rcb 
1^0.17/ V2.37/ ^5M( 


0.02 


0.4 


(19) 


Here we have fixed pb = 7 g cm“^ and 7 = 1.2 (cf. Fig¬ 
ure 3 of paper I which shows that 7 ranges from 1.2 to 
1.3 inside the innermost convective zone; although that 
figure pertains to a gas-rich nebula, similar val ues of 7 
are obtained in a gas-poor nebula). In writing (19), the 
last parameter to be calculated is the overall normaliza¬ 
tion /; the best agreement with our numerical models is 
obtained for / between 1.2 (a = 0.1 AU) and 1.3 (a = 1 
AU). (We cannot calibrate / for a > 1 AU in gas-poor 
nebulae because the relevant densities fall below those in 
our opacity tables.) 

Figures and d emo nstrate how well our semi- 

analytic scming relation ([W ) does in reproducing the full 
numerical results. We emphasize that the exponents in 


equation (19) are not merely fit parameters, but follow 
from th e physical considerations underlying equations 

(l|)-(@. 



Figure 3. Theory (dashed curves; equations |19[ [20| |22| and |24[ | 
vs. numerics (solid curves) for 5M0 cores under a variety of neo- 
ular conditions. For the most part, the agreement is good: before 
runaway, GCRs do scale with time as under many circum¬ 
stances (see also the master equation |17[ |. Solutions are truncated 
at disk depletion times: idisk,siow = lU Myr for gas-rich disks, and 
tdisk,fast = 1 Myr for gas-poor disks. The normalizations for the 
semi-analytic curves (i.e., the values of /) are adjusted by hand to 
match those of the numerical curves. 


Equation (18), o btained by numerically fitting the tabu¬ 
lated opacities of|Eerguson et al.|(|2005), defines the rele- 
vant opacity constants when evaluating equation (171 for 
dusty atmospheres. Also characterizing dusty models is 
Trcb — 2500 K: the temperature at which H 2 dissociates. 


Retaining the dependence on pj, and adopting pi, oc My.} I Va- 


lencia et al. 


2006 I yields GCR oc My/e"“A+3(i+a)/4{7-i)]/(2+^ 
I his correction hardly changes the dependence of GCR on Mcore 
for example, for dusty atmospheres, GCR oc instead of 

(see later subsections). 


Interestingly, we see in Eigurej^that the dependence 
of GCR on Mcore at fixed time t = 1 Myr differs across 
1 M 0 . Atmospheres atop core masses < IM 0 have cooled 
to their isothermal endstates and have stopped accreting 
before the sampled time; their GCRs have reached their 
maximum values. 

2.1.2. Dusty and Gas-Rich from 0.1-5 AU 

This case is almost identical to the case considered 
above. The only change in going from gas-poor (E = 
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Figure 4. GCR vs. Mcore at fixed time t = 1 Myr, demon¬ 
strating two stages of atmosphere acquisition. At low core masses 
< 1M0, GCRs have reached their maximum values: atmospheres 
have “maximally cooled” to their isothermal endstates (to cool is 
to accrete, and these planets are too cool to accrete further). Black 
crosses are numerically calculated maximum GCRs and match ex¬ 
actly the analytically computed blue curve for isothermal atmo¬ 
spheres (evaluated with T = lOOOK, the disk temperature at 0.1 
AU). At high core masses > IM0, atmospheres are still cooling and 
growing at t = 1 My r; t heir G CR s (open and filled circles) obey 
the scaling relations ( |l9[ ) and ( |20| | which predict GCR oc 
All data shown are for dusty atmospheres at 0.1 AU. 


2.3 2.8 3.6 4.8 7.5 17.0 
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o 

§ 10 “^ 
III 


0 
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z 

Figure 5. GCR vs. gas metallicity Z at fixed time t = I Myr 
for a 5M0 core in a dusty, gas-poor nebula, demonstrating that 
GCR is not a monotonic function of Z (which itself is assumed 
constant with time and space for a given model). Fo r Z < 0.2, 
the increased opacity with increased Z (equation |18^ suppresses 
cooling and yields smaller GCRs. For Z > 0.2, increases in mean 
molecular weight with Z necessitate larger Lcool to maintain hy¬ 
drostatic equilibrium; faster cooling at higher Z increases GCRs 
at a given time. The steep dependence of GCR oc predicted 

by l|19[ | leads us to compute two sets of curves: one where we 
fix Vad = 0.17 and = 2500 K (thick dashed), and another 
where we let T^cb vary according to the numerical models (thin 
dot-dashed) which naturally does better at fitting the data. 



Without the fudge factor / to mop up discrepancies, 
our derivation states that conditions at the rcb are inde¬ 
pendent of the nebular environment — i.e., T^cb = 2500 
K and Krcb = are determined by the microphysics 

governing the conversion of H 2 to H”, and p^ch as given 
by equation (11) does not depend on the outer boundary 
conditions. 'These statements are largely but not com¬ 
pletely true. In reality, there is a slight dependence of 
Prcb on the nebular density pout (as noted in paper I). 
For the same core mass, GCR, and outer boundary ra¬ 
dius, a larger pout implies a shallower atmospheric den¬ 
sity profile. Then the density at the rcb (whose tem¬ 
perature is assumed fixed at 2500 K) should be lower; 
in turn, the lower prcb reduces the optical depth and 
thereby enhances the cooling luminosity (equation 12). 
This explains qualitatively why the GCR (equivalently, 
/) is a few times larger for the gas-rich case than for the 
gas-poor case, all other factors being equal. 

Figure also illustrates how the threat of runaway gas 
accretion is greater in gas-rich disks, not only because 
they produce slightly faster cooling = slightly faster ac¬ 
creting atmospheres, but also because they last longer 
than gas-poor disks. A lOx longer lifetime enables the 
GCR to grow by an extra factor of 10°'^ = 2.5. The final 
GCR shown for the gas-rich case skirts dangerously close 
to the runaway value (formally evaluated to be 0.48; pa¬ 
per I). Although the curves for dusty atmospheres shown 
in Figure [prefer only to cores at 0.1 AU, the same 
propensity to runaway applies to dusty SMq planets at 
all orbital distances out to 5 AU (/ hardly varies between 
0.1 AU and 5 AU). _ 

Note that|Piso et al.|(|2015|) quoted a much larger crit¬ 
ical core mass of ~3OM0 at 5 AU in a dusty nebula (see 
their Figure 7). We have trac ed the origin of th e discrep¬ 
ancy to three sources. F irst, Piso et al. (2015) used the 
an alytic opacity model of Hell &: Lm (|l994|) which, unlike 
the Ferguson et al. (2005) opacities that we use, does not 
account tor dittereht sublimation temperatures of differ¬ 
ent dust species and appears to overestimate k(H“) at 
the rcb. Their higher k suppresses cooling relative to our 
models. Second, these authors defined the runaway time 
trun as the moment when Mgas/Mgas falls to 10% of its 
maximum value. This time appears systematically longer 
than our trun—defined as the time when Lcooi attains its 
minimum—by factors of 2-3. We prefer our definition as 
the minimum Lcooi has physical significance: it divides 
stable from unstable thermal equilibria in the presence 
of planetesimal acc retion (see discus sion surrounding our 
Figure 1). Finally, Piso et al. (2015) compared their t^un 
against a disk lifetime of tdisk = 3 Myr, whereas we adopt 
^disk = 10 Myr for our gas-rich models. 


Smmen/200) to gas-rich conditions (E = Emmen) is that 
the fitted normalizations are higher, running from / = 3 
(O.l AU) to 2 (1 AU) to 1.8 (5 AU): 



This scaling relation is compared against the numerical 
model in Figures]^ andthe agreement is good. 


2.2. Dust-Free Atmospheres 

Dust-free atmospheres behave qualitatively differently 
from dusty atmospheres. Removing dust as a source of 
opa city (either through grain gro wth or sedimentation; 
e.g., Mordasini|2014 OrmeH2014 ) renders the outermost 
atmospheric layers entirely radiative. The only rcb of the 
atmosphere sits at the base of this radiative and nearly 
isothermal outer shell: Trcb ~ Tout, the temperature at 
the atmosphere’s outer boundary, set by the ambient 
disk. Not surprisingly, how the GCR evolves depends 
more sensitively on nebular conditions for dust-free at- 
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mospheres than for dusty atmospheres, as the latter are 

lines are pressure broadened. At log T ~ 3, atomic alkali 

buffered by the radiative window opened by dust sub- 

metals dominate the opacity (see 

Freedman et al. 

2014 

limation (Section 2.1). Here we quote fitting formulae 

their Figure 1). 'ihe spike in n that the alkalis proc 

uce is 


for Krch and evaluate our GCR scaling relation (17) for 
dust-free atmospheres under various nebular conditions. 



log T(K) 


Figure 6. Dust-free (top) and dusty (bottom) opacities at solar 
metallicity for several densities (logp(gcm“^) = —1, -2, -10 

from top to b ottom). Black solid li nes correspond to tabulated 
opacities from |Ferguson et al.| ||2005|l, except for logp > — 6 and 
log T(K) > 3.0 where data are extrapolated following the method 
described in paper I. Blue dashed li nes correspond to tabulated 
opacities from|Freedman et al.|||2014|l, except for logp = —10 and 
log T <2.5 where tt i eir an alytic tits are used (their equations 3- 
5). |Freedman et al.| (|2014|l present their opacity tables in (P, T) 
space so we use p = 2.;3Y4 to convert P to p. The abrupt jumps in 
the dust-free k at T = 1000 K for logp ^ — 3 are due to a failure 
in equation-of-state calculations (J. Ferguson, private communica¬ 
tion). Our model atmospheres never attain such high densities at 
T = 1000 K so this bug is irrelevant. 

The opacities under lying our models are drawn from 
Ferguson et al.| (|2005 1 and Freedman et al. (2014). We 
employ the forrner for log > 2.7 (the domain of 

their tabulations), and the latter for colder temperatures. 
A smooth merging of the tw o datasets is effected by off- 


Freedman et al. 


(2014) to 
005p data at log V = 2.7. 


setting all of the data from 
match the |Ferguson et al.| 

Depending on p, the otlset in k ranges from 0.1 to 0.4 
dex. If no numerically tabulated value for k is available 
for a desired (p, T) below log T = 2.7, we use the ana¬ 
lytic f its given by equations (3)-(5) of Freedman et al. 
(2014); above log T = 3.6, we extrapolate following the 


procedure described in paper 1 . 

The top panel of Figure plots the dust-free k{p,T). 
We surmise the following features. At log T < 3, the 
main absorbers are water, methane, and ammonia whose 


washed out at large p, presumably from pressure broad¬ 
ening. At still higher log T > 3.3, H“ reigns, as it does 
for dusty atmospheres. 

Two-dimensional power laws fitted to k{p,T) in re¬ 
stricted domains of (p, T) are p resented below. I t is as- 
sumed throughout that k oc Z; |Freedman et al.| (|2014| 
found a nearly linear Z-dependence tor 1' ~ 25U-30UU K. 

2.2.1. Dust-Free and Gas-Poor at 0.1 AU 
For dust-free atmospheres at 0.1 AU, 

«:(dust-free, 0.1 AU) ~ 7 x 10“^ cm^ g“^ 
p \ ( 1 \ f Z 


10“® g cm“3 


1000 K 


0.02 ’ 


valid for —5 < log p(g cm ^) < —3 and 1000 < T(K) < 
2500. The corresponding GCR is 


GCR ~ 0.16 ( /- 

-1.. O 


0.02 


t 


0.4 


1 Myr 

Vad 

0T7 


0.4 


1600 K 


V ^rcb 
/ Mrcb \ 

v^y 


1.9 


3.3 fM, 


5ilTr 


1.6 


( 22 ) 

where the nominal values for Vad and T^ch are drawn 
from our full numerical model at 0.1 AU (for comparison, 
Tout = 1000 K), and where we have calibrated / = 1.3. 
Equation ( [2^ is plotted against the full numerical solu¬ 
tion in Figure 

The K quoted above is most relevant for gas-poor disks 
at ~0.1 AU whose surface densities E = Emmen/200. 
The opacity behaves differently in gas-rich disks. We 
do not present results for dust-free atmospheres in gas- 
rich disks at 0.1 AU because such atmospheres stay fully 
convective for GCRs at least up to ~0.3 and cannot be 
evolved using our numerical model. There would not be 
much point to following them anyway, since the combi¬ 
nation of dust-free and gas-rich conditions leads to rapid 
runaway. 

2.2.2. Dust-Free and Gas-Rich Beyond 1 AU 

Outside 1 AU, temperatures fall. For 100 < r(K) < 
800 and —6 < log p(g cm“^) < —4, 

K(dust-free, > 1 AU) ~ 1 x 10“® cm^ g“^ 

^ \ 0.6 / p \ / Z ^ ^ 

10-6g cm-3 ) I^IOOkJ 


(23) 


whence 

GCR^ 


0.1 (T 
2.8 


t /2ooKy® 

Ikyry V Trcb J 


/y^y" /^y-2 / Meou 
Z J \O.25) V2.37/ 


(24) 


where the nominal Vad = 0.25 (equivalently, 7 = 1.3) 
represents a rough average over our dust-free numerical 
models at distances > 1 AU. The numerical models also 













































indicate that at these large distances, Trcb nearly equals 
Tout; more distant plan ets have more nearly isothermal 
outer layers. Equation (241, normalized for nebular con¬ 
ditions at 5 AU, is testea against the full numerical so¬ 
lution in Fi gure [3| 

Equation (24) highlights just how susceptible dust-free 


super-Earths are to runaway at large distances in gas-rich 
disks. At 5 AU, T^ch = Tout — 200 K and a core will 
attain GCR ~ 0.5 in as short a time as 0.05 Myr. (We do 
not present solutions for gas-poor disks at large distances 
because the relevant gas densities fall below those covered 
by our opacity tables. It is easy, however, to guess what 
these solutions would look like: GCRs would increase 
on timescales almost as fast as those in gas-rich disks, 
and would be limited only by the gas available and the 
remaining disk lifetime.) 

3. SUMMARY 


We have developed an analytic model describing how 
a solid core accretes gas from its parent disk. The model 
assumes that the planet’s growing atmosphere cools pas¬ 
sively: heat inputs from external sources such as plan- 
etesimal accretion are ignored. Heating from planetesi- 
mal accretion can be relevant for low-mass < 5 M 0 cores, 
but only under exceptionally fine-tuned conditions. Too 
high a planetesimal accretion rate leads to excessively 
massive cores which undergo runaway, while too low an 
accretion rate is energetically irrelevant. The window 
for planetesimal accretion to be relevant is so narrow — 
it vanishes completely for core masses > 5M(^, which 
cannot help but undergo runaway in gas-rich disks that 
persist for ~10 Myr — that neglecting energy deposition 
by solids is almost certainly safe. 

We derived how the gas-to-core mass ratio (GGR) 
varies with time t, core mass Mcore, and metallicity Z 
for passively cooling at mos pheres. The scaling relation¬ 
ship, given by equation (ITtI), applies in various formation 
environments at times preceding the onset of runaway 
gas accretion (at times when the envelope’s self-gravity 
is still negligible, i.e., when GCR < 0.5). The analytic re¬ 
sults agree with numerical calculations to within factors 
of 1-3. 

Three physical ingredients make our derivation possi¬ 
ble. (1) To cool is to accrete: in the absence of any heat 
source, atmospheres of mass Mgas grow on a timescale 

equal to their cooling time: Mgas/Mgas ~ \E\/Lcooi ~ t, 
where E is the atmosphere’s energy and Lcooi is its cool¬ 
ing luminosity. The final GCR is then set by the condi¬ 
tion that the atmosphere’s cooling time tcooi equals the 
gas disk’s depletion time tdisk- (2) Most of the mass of 
the atmosphere is concentrated toward the core, whose 
given properties provide an easy reference for scaling our 
variables. The mass concentration follows from density 
gradients made steep by the dissociation of H 2 and the 
consequent lowering of the adiabatic index 7 down to 
1.2-1.3. (3) The cooling luminosity is regulated by the 
radiative-convective boundary (rcb) and therefore can be 
evaluated there. Specifying the microphysical properties 
of this boundary (temperature, opacity) enables us to 
quantify the rate of cooling. 

The scaling indices of GCR (t, Mcore, Z) are deter¬ 
mined by the innermost convective zone’s adiabatic in¬ 
dex 7 (calibrated using our numerical models) and the 
dependencies of the opacity at the rcb. Dusty atmo¬ 


spheres have rcb’s that are distinct from those of atmo¬ 
spheres rendered dust-free by grain gro wth and sedimen¬ 
tation. In dusty atmospheres (Section 2.1), the rcb oc¬ 
curs universally — at all orbital distances and in gas-rich 
or gas-poor nebulae — at the H 2 dissociation front, set¬ 
ting the temperature T^cb — 2500 K. The rcb in dusty 
atmospheres is insensitive to external nebular conditions 
because dust sublimation opens a radiative window in 
the interior of the planet, whose base (the rcb) is deter¬ 
mined by the local microphysics of how H 2 converts into 
H“ whose opacity rises sharply with temperature. All 
dusty atmospheres evolve as GCR oc 
the scaling indices apply globally. Note that the mean 
molecular weight /i depends inversely on (1 — Z) so that 
all other factors being equal, GCR first decreases with Z 
(as Z~^ '^; this dependence arises from opacity) and then 
increases (via /i^'^) as atmosphere s bec ome “heavier”. 

Dust-free atmospheres (Section 2.2) are more trans¬ 
parent in their outer layers and their rcb’s are lo¬ 
cated at higher altitudes: conveniently, Tjjutj 

the given temperature of the ambient nebula. For 
dust-free atmospheres, the dependence of rcb proper¬ 
ties on nebular properties implies that the GCR scal¬ 
ing indices change with orbital distance, depending on 
the behavior of the local opacity. We examined two 
cases. At 0.1 AU in a gas-poor nebula, GCR oc 


Beyond 1 AU in a gas-rich neb¬ 


ula, GCR cx 


cases 

we evaluated, GCR oc The GCR-t scalin g is de¬ 

termined by the Krcb-Prcb Scaling (see equation 17). It 
appears the latter scaling does not change much across 
environments (/trcb oc where a ~ 0.3-0. 6 ; this ap¬ 
pears to be an average between a = 1 characterizing 
pressure-broadened opacities at T ^ 100 K and a « 0 
ch aracterizing the opaci ty at T ^ 2000 K; see Figure 7 
of Freedman et al.|[2014 ). 


While dusty atmospheres behave more-or-less the same 
way under a variety of nebular conditions, dust-free at¬ 
mospheres depend more sensitively on disk temperatures. 
Gas opacities tend to decrease with colder temperatures 
and therefore dust-free atmospheres grow faster the far¬ 
ther they are from their central stars. 

Atmospheres accrete faster in gas-rich environments 
than in gas-poor ones, but not by much. For example, 
in our numerical models of dusty atmospheres, dropping 
the nebular density by a factor of 200 drops the GCR 
by a factor of 2.5, all other factors being equal. We 
have checked that similar results obtain for dust-free at¬ 
mospheres. Without introducing a no rmalization correc¬ 
tion (our factor of / in equation [l7| ) , our analytic scal¬ 
ings for dusty atmospheres would predict no depen dence 


at all on the ambient nebular density (see Section 2.1.2 
for the technical reason why in reality there is a small 
dependence). The insensitivity to nebular density sug¬ 
gests that the accretion of planetary atmospheres pro¬ 
ceeds about the same way whether or not planets open 
gaps in disks. Gap opening in viscous gas disks does 
not starve the planet; it is well-recognized that material 
continues to flow past the planet, through th e gap (e.g., 
Lubow fc D’Angelo|2006 Duffell et al.|2014 ). Moreover, 


the factors by which densities are suppressed in gaps are 


likely overe stimated in 2D simulations (e.g., Kanagawa 


et al. 2015): 3D gaps are messier. Careful studies of 
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the 3D flow dynamics of planets embedded in disks (e.g., 


D^mg elo fc Bodenheimer|2013[ |Ormel et al.||2015 Fung 

et a] '.|2015f are clearly needed to test whether gap open- 
ing and/or hydrodynamic effects substantively alter the 
ID accretion theory we have presented. 

The analytic solutions for gas-to-core ratios developed 
here enable us to calculate the properties of planetary 
envelopes just before protoplanetary disk gas disperses. 
The solutions — encapsulated in the condition tcooi ~ 
idisk — provide initial conditions for models of subse¬ 
quent atmospheric loss by hydrodyna mic winds powered 
by st ellar irradiation or internal heat (Owen & Wu|2013 
2015 1. Our model complements studies that include de¬ 
tailed molecular ch emistry (e.g., Hori & Ikoma 2011 Ven- 
turini et al. 2015); we find good agreement with those 
modefs tor 'Z up to ~0.5. The simplicity of our solutions 
make them suitable for the construction and/or diagnosis 
of more rea listic models that include 3D hydrodynamic 
effects (e.g., D’Angelo fc Bodenheimer|2013 Ormel et al. 
2015 Fung et al. 2015). I'he solutions can also be used 
to quickly assess the plausibility of various planet forma¬ 
tion scenarios. This is the task we set ourselves for paper 
III. 
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